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We solve numerically equations of motion for real self-interacting scalar fields in the background of 
Reissner-Nordstrom black hole and obtained a sequence of static axisymmetric solutions representing 
thick domain walls charged black hole systems. In the case of extremal Reissner-Nordstrom black 
hole solution we find that there is a parameter depending on the black hole mass and the width of 
the domain wall which constitutes the upper limit for the expulsion to occur. 

I. INTRODUCTION 

Relics of cosmological phase transitions are considered to play an important role in cosmology [1]. They are signs 
of the high-energy phenomena which are beyond the range of contemporary accelerators. Such topological defects as 
cosmic strings attracted great interests due to the motivation of the uniqueness theorems for black holes and Wheeler's 
conjecture black holes have no hair. Black hole cosmic string configurations were studied from different points of view 
[2] revealing that cosmic string could thread it or could be expelled from black hole. As far as the dilaton gravity was 
concerned it was also justified [3,4] that a vortex could be treated by the remote observer as a hair on the black hole. 
On the other hand, for the extreme dilaton black hole one has always the expulsion of the Higgs field (the so-called 
Meissner effect) . The problem of vortices in de Sitter (dS) background was analyzed in Ref . [5] , while the behaviour 
of Abelian Higgs vortex solutions in Schwarzschild anti de Sitter (AdS), Kerr, Kerr-AdS and Reissner-Nordstrom AdS 
was studied in Refs. [6,7]. The vortex solution for Abelian Higgs field Eqs. in the background of four-dimensional 
black string was obtained in Ref. [8]. 

In order to investigate the behaviour of the domain walls or cosmic strings configurations in curved spacetimc a 
number of works use a thin wall or string approximations, i.e., infinitely thin non-gravitating membrane described 
by the Nambu-Goto Lagrangian. Dynamics of scattering and capture process of an infinitely thin cosmic string in 
the spacetime Schwarzshild black hole were elaborated in Ref. [9]. Christensen et al. [10] revealed that there existed 
a family of infinitely thin walls intersecting black hole event horizon. Then, the considerations were applied to the 
scattering of inifitely long cosmic string by a rotating black hole [11]. 

Nowadays, the idea that the Universe is embedded in higher-dimensional spacetime acquires much attention. The 
resurgence of the idea is motivated by the possibility of resolving the hierarchy problem [18], i.e., the difference in 
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magnitudes of the Planck scale and the electroweak scales. Also the most promising candidate for a unified theory 
of Nature, superstring theory predicts the existence of the so-called D-branes which in turn renews the idea of brane 
worlds. From this theory point of view the model of brane world stems from the idea presented by Horava and 
Witten [19]. Namely the strong coupling limit of E$ x E 8 heterotic string theory at low energy is described by 
eleven-dimensional supergravity with eleventh dimension compactified on an orbifold with Z 2 symmetry. The two 
boundaries of spacetime are ten-dimensional planes, to which gauge theories are confined. Next in Refs. [20,21] it 
was argued that six of the eleven dimensions could be consistently compactified and in the limit spacetime looked 
five-dimensional with four-dimensional boundary brane. 

Then, it is intriguing to study the interplay between black holes and domain walls (branes) . Nowadays, this challenge 
acquires more attention. The problem of stability of a Nambu-Goto membrane in Reissner-Nordsrom de Sitter 
(RN-dS) spacetime was studied in Rcf. [12], while the gravitationally interacting system of a thick domain wall and 
Schwarschild black hole was considered in [13,14]. Emparan et al. [15] elaborated the problem of a black hole on 
a topological domain wall. In [16] the dilaton black hole-domain wall system was studied analytically and it was 
revealed that for the extremal dilaton black hole one had to do with the expulsion of the scalar field from the black 
hole (the so-called Meissner effect). The numerical studies of domain wall in the spacetime of dilaton black hole was 
considered in Ref. [17], where the thickness of the domain wall and a potential of the scalar field ip 4 and sine-Gordon 
were taken into account. 

In our paper we try to provide some continuity with the previous works [16,17] concerning the behaviour of dilaton 
black hole domain wall systems, namely we shall consider the thick domain wall in the background of RN black hole. 

The paper is organized as follows. Sec. II is devoted to the basic equations of the considered problem. In Sec. Ill we 
studied the RN-domain wall system both for cp 4 and sine-Gordon potential. We pay special attention to the extremal 
RN black hole case and reveal that there is an expulsion of the domain wall fields under certain conditions. In Sec. IV 
we concluded our investigations. 

II. REISSNER-NORDSTROM BLACK HOLE DOMAIN WALL SYSTEM 

In our paper we shall consider the spherically symmetric solution of Einstein-Maxwell equations described by the 
RN black hole metric of the form 

ds 2 = - ( 1 - — + 91) d t> + + r 2 {d6 2 + sin 2 0d<p 2 ). (1) 



The metric is asymptotically flat in the sense that the spacetime contains a data set (S e nd, gij,Kij) with gauge fields 
such that T, end is diffeomorphic to R 3 minus a ball and the following asymptotic conditions are fulfilled: 



\9H ~ h\ +r\da9ij\ + ■■■ + r k \d ai ... ak g ij \+r\K ij \ + ... + r k \d ai ... ak K tJ \ < of^J, (2) 

\F afj \+r\d a F afj \ + ... + r k \d ai ... ak F a0 \<O^. (3) 

RN black hole spacetime will be the background metric on which we shall consider the domain wall equations of 
motion built of a self-interacting scalar field. In our investigations we use a general matter Lagrangian with real Higgs 
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field and the symmetry breaking potential of the form as follows: 



jC dw = --V^Wip-U(ip), (4) 



where ip real Higgs field and the symmetry breaking potential U(ip) has a discrete set of degenerate minima. The 
energy-momentum tensor for scalar field may be written as follows: 

Tij(<p) = ~^9ij V m yA7"V -U(ip)gij + VitpVjip. (5) 

It is turned out that for the convenience we can scale out parameters via transformation X = ip/rj and e = S-kGtj 1 
[15]. The parameter e is responsible for the gravitational strength and it is also connected with the gravitational 
interaction of the considered Higgs field. One can also define the potential V(X) = ^y^-, where Vp = Ar/ 4 . All these 
quantities enable us to rewrite the Lagrangian for the domain wall fields in a more suitable form. Namely it now 
reads 



SirGCdw = k 



w2 ^xwx + v{x) 



(6) 



where w = Sv q Vf represents the inverse mass of the scalar after symmetry breaking. It characterizes the width of 
the wall defect. Just using Eq.(6) we arrive at the following expression for X field: 

1 dV 

v ^ x -^ox-°- W 

As in Rcf. [17] we shall take into account for two cases of potentials with a discrete set of degenerate minima, namely 
the ip 4 potential described by the relation 

U^) = \{^-rf)\ (8) 

and the sine-Gordon potential in the form of 

U 2 (<p) = Ar/ 4 [1 + cos^/77)] . (9) 

III. BOUNDARY CONDITIONS 
A. Boundary conditions for <p 4 potential 

As we mentioned RN black hole spacetime is asymptotically flat, thus the asymptotic boundary solution of equation 
of motion for potential (8) will be the solution of equation of motion in flat spacetime. It reads as 

ipxiz) = r/tanh( v /A72?7^), (10) 
or rewritten in our units implies the following: 

X M )=tanh(^£) (11) 
Thus, the equation of motion (7) for the scalar field X in the considered background has the form of 
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and the energy density yields 
E 



2Mr + Q 2 )d r X 







+ 2 • O de 


sin QdeX 


r z sin 9 





-— X{X 2 -1) = 0. 



w 2 -i(A 2 -l) 2 . 
4 



2 V 7 V r r 2 7 2 
On the horizon of the RN black hole relation (12) gives the following boundary condition: 



■a. A" I 



1 



[sm9d e X] + ^X{X 2 -1) 



(12) 



(13) 



(14) 



lr=r + r^sinfl 

where r + = M + \J M 2 — Q 2 and r_ = M — a/M 2 — Q 2 are respectively the outer and the inner horizon of the RN 
black hole. 

As in Ref. [17] we shall investigate only the case when the core of the wall is located in the equatorial plane 9 = n/2 
of the black hole, we impose the Dirichlet boundary condition at the equatorial plane 



(15) 



as well as the regularity of the scalar field on the symmetric axis requires the Neumann boundary condition on the 
z-axis: 



dX 



0. 



(16) 



8=0 



For the asymptotical flatness of the RN black hole solution we should have far from the black hole the flat spacetime 
solution (10). But we are limited by our computational grid which is finite, then we require the following boundary 
condition ought to be satisfy at the boundary of the grid 



tanh 



^m.ax COS ( 



W 



(17) 



B. Boundary conditions for the sine-Gordon potential. 



As far as the sine-Gordon potential (9) is concerned the flat spacetime solution is given by 



^2(2) = r\ |4 arctan 
or in our units this is equivalent to following: 

X(r, 9) = 4 arctan 



cxp (y/\i]z) — 7r| 



exp 



In this case V(X) = 1 + cos(A) and dV/dX = — sin(A), and thus the equation of motion (7) takes the form 



^ d r (r 2 - 2Mr + Q 2 J d r X 
while the energy density is given by the relation 



+ 2 ■ a de 
sin 9 



sin 9de X 



(18) 



(19) 



(20) 
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E = 



-(d r X) 



2M Q 
+ — r 



1 - cos(X). 



On the horizon from Eq.(20) one has the boundary condition as follows: 



-d r X\ 



1 



r\ sin ( 



[sinddeX] ^ sin(X). 



(21) 



(22) 



Of course, this must be accompanied by the Dirichlct boundary conditions at the equatorial plane of the black hole 
(15), and the Neumann boundary condition on the symmetry axis (16), and at the outer edge of the grid 

r max COS 9 



X 



= 4 arctan 



exp 



w 



— 7T. 



(23) 



Another interesting problem that should attract attention to is the case of extremal RN black hole domain wall 
configuration. As was revealed analitycally in Rcf. [15] that for some range of parameters we should have the 
expulsion of the domain wall, contrary to the situation which takes places in extremal dilaton black holes [16,17]. 
Simple arguments given in Ref. [15] it was revealed that for at least very small extremal black hole sitting inside the 
domain wall, the black hole would expel it. The mass bound for the extremal RN black hole for the expulsion was 
also found. For w = 1 the black hole mass have to satisfy the condition M 2 < \. 

When r + = r_ we have to do with the extremal RN black hole and it implies the condition Q — M. For this case 
the metric has the following form: 

ds 2 =-(^l-^j 2 dt 2 + - ^ 2 + r 2 (d9 2 + sin 2 9d<p 2 ). (24) 

("^ r ) 

Equation of motion for the domain wall's scalar fields in the background of RN extremal black hole yields 



(r-M) 2 d r X 







+ 2 ■ a de 


sin Ode X 


r z sm 





while the relation for energy implies the following: 



E 



-Sdrxyii-- 



M 



w 



1 dV 
^OX' 



V(X). 



(25) 



(26) 



IV. NUMERICAL CALCULATIONS 
A. Reissner-Nordstrom black hole 

The method of solving equations of motion for the scalar fields X is the same as in Ref. [17]. We applied the 
modified overrelaxation method in order to take into account the boundary conditions on the considered black hole 
horizon. Eqs. of motion will be solved on a uniformly spaced polar grid (r,,^) with the boundaries at r min = r + , 
the outer radius will be r max 3> r + (we take usually r max = 20r+ ). The angle 9 changes from to n/2. Then, 
one approximates the derivatives with finite-difference expressions on the grid. The rest of the solution for scalar 
field X will be obtained from the symmetry condition. The example solutions of the equation of motion for different 
parameters and field configurations are shown on Figs 1, 2, 3, and 4. 
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B. An extreme Reissner-Nordstrom black hole 



For the extremal RN black hole scalar fields on the horizon decouple from the rest of the grid and equation of 
motion (25) becomes an ordinary differential equation: 

1 r 2 BV 

doeX + — d e X - -± = 0. (27) 
tant> w z aX 

For the reasons given in [2,17] one must carefully examine the solution on the horizon. The solution X = 0, which 
corresponds to the expulsion of the field from the black hole (the Meissner effect) always solves relation (27). But 
sometimes other solution is possible and energetically more favourable. The previous studies envisage ( [2,17]) that 
the relative size of the domain wall compared to the size of the black hole becomes a factor which determines the 

r 2 

overall behaviour of the field on the horizon. For this reason we parameterize equation (27) with parameter a = 
and solve (27) numerically using two point boundary relaxation method ( [22]). Fig. 5 shows the dependence of the 
scalar field X on the black hole z-axis as a function of the parameter a for the (p* potential. Fig. 6 depicts the same 
dependence for the sine-Gordon potential. The difference of the scalar field behaviours is clearly visible. For the value 
a < 2.0 we have the field expulsion, while for a > 2.0 the field pierces the horizon of the extremal RN black hole. 
The value of the critical parameter is the same in both cases, which is easily understood when one considers the series 
expansion of the dV/dX around zero and neglects higher terms. 

Figs 7, 8, 9, 10 show examples of the field configuration for various combinations of the a parameter and the 
adequate potentials. 

Summing it all up we conclude that the numerical studies of the scalar domain wall field in the vicinity of RN black 
hole reveal that for the extremal case one has not always expulsion of the fields from the black hole. This phenomenon 
depends on the value of the black hole radius r + and the width of the domain wall. If the parameter a < 2.0 we have 
expulsion of the scalar field X while for a > 2.0 the scalar fields penetrate the extremal RN black hole horizon. The 
domain wall extremal RN black hole system was studied analytically in Rcf. [15]. By simple argument it was revealed 
that at least small extremal RN black holes sitting inside the domain wall would expel it. Then a more precise bound 
on the mass of the considered black hole was established. For M 2 < 1/2 the defect's field cease to penetrate the black 
hole horizon, but for M 2 > 1/2 one has no such expelling. In Rcf. [15] the considerations were conducted for the case 
of the domain wall's width equaled to 1. 

On the contrary, in the case of extremal dilaton black hole domain wall system it was shown that one had always 
expulsion of the scalar field [16,17]. It turned out that the simplest generalization of Einstein-Maxwell theory by 
adding massless dilaton field dramatically changed the structure of the extremal black hole. The similar situation 
takes place when one consideres anothe topological defect, i.e., cosmic string. For extremal dilaton black holes we 
observe always expulsion of the fields [3] . 

V. CONCLUSIONS 

In our paper we elaborated the problem of thick domain walls in the background of a spherically symmetric 
solution of Einstein-Maxwell field equations. The domain wall equations of motion were built of a self-interacting real 
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scalar fields with the symmetry breaking potential having a discrete set of degenerate minima. Using the modified 
overrelaxation method modified in order to comprise the boundary conditions of the scalar field X on the black hole 
horizon we solved numerically Eqs. of motion for ip 4 and sine-Gordon potentials. All our solutions depended on the 
parameter w — 1/VXrj which was responsible for the domain wall thickness. We revealed the systems of axisymmetric 
scalar field configurations representing thick domain wall in the nearby of RN black hole for both kinds of potentials. 
We paid special attention to the extremal RN black hole solution. We parametrized the equations of motion with the 

r 2 

parameter a — and solved it by means of the two-point boundary relaxation method. For a < 2 one observed the 
complete expulsion of the scalar field X from the extremal RN black hole (the so-called Meissner effect), while for 
a > 2 the domain wall's field penetrated the extremal RN black hole horizon. This behaviour is in contrast to the 
behaviour of scalar field in the background of the extremal dilaton black hole, where we have a complete expulsion of 
the scalar fields. 
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X X 



FIG. 1. The field X (left panels) and the energy E (right panels) for the 4> 4 potential and the Reissner-Nordstrom black hole. 
Isolines on bottom panels are drawn for 0.2, 0.4, 0.6 and 0.8 for the field X and for —0.1, —0.2, —0.3 and —0.4 for the energy. 
Inlets in bottom plots show the value of the fields on the black hole horizon. Black hole has M = 1.0, Q = 0.1 and the domain 
width is w — 1.0. 
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FIG. 2. The field X (left panels) and the energy E (right panels) for the 4> 4 potential and the Reissner-Nordstrom black hole. 
Isolines on bottom panels are drawn for 0.2, 0.4, 0.6 and 0.8 for the field X and for —0.1, —0.2, —0.3 and —0.4 for the energy. 
Black hole has M = 1.0, Q = 0.1 and the domain width is w = 10. 
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FIG. 3. The field X (left panels) and the energy E (right panels) for the sine-Gordon potential and the Reissner-Nordstrom 
black hole. Isolines on bottom panels are drawn for 0.2n, OAn, 0.6tt and 0.8-7T for the field X and for —0.5, —1.5, —2.5 and —3.5 
for the energy. Inlets in bottom plots show the value of the fields on the black hole horizon. Black hole has M = 1.0, Q — 0.1 
and the domain width is w = 1.0. 
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FIG. 4. The field X (left panels) and the energy E (right panels) for the sine-Gordon potential and the Reissner-Nordstrom 
black hole. Isolines on bottom panels are drawn for 0.2n, OA-k, 0.6-7T and 0.8-7T for the field X and for —0.5, —1.5, —2.5 and —3.5 
for the energy. Black hole has M = 1.0, Q = 0.1 and the domain width is w = 10. 
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FIG. 5. Value of the field X on the black hole axis (8 = 0) as a function of the parameter a for extreme Reissner-Nordstrom 
black hole and <f> 4 potential. 
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FIG. 6. Same as Fig. 5, but for the sine-Gordon potential. 



14 




X X 



FIG. 7. The field X (left panels) and the energy E (right panels) for the (j) 4 potential and the extreme Reissner-Nordstrom 
black hole. Isolines on bottom panels are drawn for 0.2, 0.4, 0.6 and 0.8 for the field X and for —0.1, —0.2, —0.3 and —0.4 for 
the energy. Inlets in bottom plots show the value of the fields on the black hole horizon. Black hole has M = 1.0, Q = 1.0 and 
the domain width is w = 1.0, so the parameter a = 1.0. 
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FIG. 8. The field X (left panels) and the energy E (right panels) for the (j) 4 potential and the extreme Reissner-Nordstrom 
black hole. Isolines on bottom panels are drawn for 0.2, 0.4, 0.6 and 0.8 for the field X and for —0.1, —0.2, —0.3 and —0.4 for 
the energy. Inlets in bottom plots show the value of the fields on the black hole horizon. Black hole has M = 1.0, Q = 1.0 and 
the domain width is w = V0A, so the parameter a = 10.0. 
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FIG. 9. The field X (left panels) and the energy E (right panels) for the sine-Gordon potential and the extreme Reiss- 
ner-Nordstrom black hole. Isolines on bottom panels are drawn for 0.2n, OAn, 0.6n and 0.8-7T for the field X and for —0.5, —1.5, 
—2.5 and —3.5 for the energy. Inlets in bottom plots show the value of the fields on the black hole horizon. Black hole has 
M = 1.0, Q = 1.0 and the domain width is w = 1.0, so the parameter a — 1.0. 
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FIG. 10. The field X (left panels) and the energy E (right panels) for the sine-Gordon potential and the extreme Reiss- 
ner-Nordstrom black hole. Isolines on bottom panels are drawn for 0.2n, OAn, 0.6n and 0.8-7T for the field X and for —0.5, —1.5, 
— 2.5 and —3.5 for the energy. Inlets in bottom plots show the value of the fields on the black hole horizon. Black hole has 
M = 1.0, Q = 1.0 and the domain width is w = \/6TT, so the parameter a = 10.0. 
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